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The following list of changes, which includes the changes listed in 
all previously issued errata, should he no teds 


Page 3, line 5» Read 0 instead of A, 

Page 7, equation (6). The expression in the parentheses, last line, 
should read; 


( 


m V 1-at 2 


EMI 4. 


0 


read: 


Page 8, equation (8). The expression in brackets, last line, should 


mtVl-a-t 2 ~ (mt-at) 


Page 13. The definition of £ should read: 


t - - 3 £- 

1+S 2 € a 


Page 14, line preceding equation (19) » 
to read: 


Correct expression for x a 


x a 



Page IT, line preceding equation (28), Read "wing" for "tip.” 

Page 18, equations (31) and (32). The right-hand sides should be 
multiplied by 2 (i.e., the number 2 appearing in the coefficients should 
be 4). 
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Page 20, equation ( 36 ) . Insert p in the denominator of the 
coefficient of the right-hand member. 

Page 20, bottom of page. The line preceding the series for K(k* ) 
should read , . when k is small.” 

Page 23/ equations (46) and (4-7) . Multiply right-hand members 

■by (-1). 

Page 24, equation for H' (m) . Change to read? 


m-a t 

a t BHe t J 

r x - ( 

l\ f lo-ax. ' 

\ 

a t a m 4 (l+m) 

m 2a-t(l+m) ' 

L V 

» 2m/ L * 6( l+m) - 

j 


Page 26, equation (51). Replace by 


u = r.p 


F(<p, «/l-^t 2 ) 

K( 


where K and P are the co mplete and incomplete elliptic integrals of 
the first kind, of modulus ,/l— m^ 2 , and 

(p * sin-i / 1 ~ l o . 

V 1-Ha t 2 

(This solution is due to H. S. Ribner, of the NACA.) 

Page 27, equation (55). The expression in brackets should read: 

\ _ */ g 

- " K( 

p 

Page 28, lines 2, 3, 4, and 5. The series expansion of 1 E is 

1 n 

no longer pertinent. Substitute the following for lines 2, 3/ 4, and 5* 


The following series expansion (from reference 9) is useful; 

x ffi— -fx.— 1— . . 

KWl-“t 2 ) SK( L V X 2 % 2 6 2 


where 
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Page 28, equation (57)- Change to read: 


AL 


— if-m^s 2 


qa a t 2 (m t 2 -® 2 )E( */l-m 2 ) 


1_ n^rn^.— a^) 


(maHn t a t ) Vm^t 2 I " 1 - - 

L v m t“ a t. 

^m t (l-^i) K+m) || j^(m+a t ) 2 + (nh-a t ) 2 ^ 


[ F($ t ,k) - 


} 


F($oA)] 2 —f [(m 2 -m t a t ) 2 +3m 2 (m t -a t ) 2 ] [E(4>t,k)-E($ 0 ,k) ] 

HLj- —iH / 


Page 29> equation ( 58 ). Close trace at end of equation. 

Page 30. The expression for the induced moment AM should read; 


AM 


- lt|> y |3B,1 ° r 1 ( . J-\ F (t, vr^5) i®. at< 
K( v^at 5 ) 'An t V at m t 3to/ dt 0 


and 


AM _ 4m3s3 


r “f^t 

r, ”/a 1.1 

", B( /l~^t 2 ) 

) 

l a t 

L K( ^l-m t 2)J 3 

- K( /y/l-rn,- 2 )- 

j 


( 60 ) 


Page 3^, line 21. Amend to read "(Au) a replaced by the right-hand 
member of equation ( 51 )” and delete the expression -which follows. 

Page 38 * following equation (74) . Read "equation ( 38 ) " for 
"equation ( 30) . " 
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Page 40, equation (77). The expression in "brackets should read; 

h§ 

m t K( v 1-04.2) 

Page 46. The table should read as follows: 


Wing 

Aspect 

ratio 

Glee 

Cd/Cl 2 

Tapered wing 

3.85 

3.09 

0.160 

Untapered wing 

1.73 

1.-97 

.326 

Modified untapered 
wing 

2.45 

2.62 

.236 


Page 49. In e, read p 2 z 2 for pz 2 

In g(a), read for & 
m t m 

In Q, read (1— n^) for ( l+m t ) 

Page 52, equation (A6). Insert a 4- 2 in denominator of the coeffi- 
cient of the right-hand member. 

Figures 

The changes arising from the correction of equation (51) have only 
a small effect on the calculated results presented. (There is a major 
change in the downwash, which was not evaluated in this report.) The 
general shape of the pressure distribution remains the same. The 
magnitude of the correction is indicated in the figure below, in which 
the correct Au/u 0 and the curve given by the original equation (51) 



-1.0 -m t 0 m, iO 


to 




NACA TN 1555 


ERRATA 


5 


are compared for m^=0.57- This figure corrects figure 13, and may he 
used to find the percent error in component 1 of the pressure distribu- 
tions in figure 17 . Component 5 in figure 17(d) will likewise be 
increased. 


Tables I and II 

Table I. The component of — calculated for the untapered wing 
from equation (55) should be changed from — 15-90 to —18.28, and from 
equation ( 73 ), from +0.14 to +0.16. The total is therefore 114.28 
and 01^=1.97, or 0.034 per degree. 

M 

Table II. The values of g-j contributed by the symmetrical wake 
solution, as corrected, are as follows: 19-9/ 207. 3, and 4l.3« The 

amended totals are —457.4, —728.2, and — 66o.5» and the ratios of moment 
arms are 1 . 08 , 0 . 82 , and 1.01. 


NACA-Langley - 8-8-49 - 1000 
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TECHNICAL NOTE NO. 1555 

THE THEORETICAL LIFT OF FLAT SWEPT-BACK 
WINGS AT SUPERSONIC SPEEDS 
By Doris Cohen 


SUMMARY 

The method of superposition of linearized conical flows has 
"been applied to the calculation of the aerodynamic properties of 
thin, flat, swept— hack wings at an angle of attack. 

Various cases are distinguished, depending on the sweep of the 
leading and trailing edges relative to their respective Mach lines. 
Where the Mach line from the apex of the wing lies behind the 
leading edge hut intersects the tip, formulas for the total lift 
are given in closed form. The induced drag is simply the lift times 
the angle of attack. To obtain the pitching moment the numerical 
integration of elementary functions is required. Where the leading 
edge, hut not the trailing edge* is swept behind the Mach lines, 
the pressure distribution and the total lift are both given in 
closed form, but not the pitching moment. The induced drag is 
calculated by a simple formula from the lift. 

The highly swept wing with both leading and trailing edges 
behind their respective Mach lines cannot be completely solved by 
the present method except by a series of successive approximations. 
However, the pressure distribution can be determined over all but 
the generally small region in which the tip and trailing-edge Mach 
cones overlap, and the total lift, pitching moment and induced drag 
can be approximated with practical accuracy, provided the aspect 
ratio is not so high nor the Mach number so low that the Mach lines 
from the trailing edge intersect the leading edge. (The latter case 
is not treated in the present paper.) 

The complete procedure for the highly swept wing is outlined 
and various degrees of approximation suggested. Application of 
the general method to other plan forms is also discussed. 

The investigation has produced an interesting result regarding 
the lift near the tip of a highly swept wing. The formulas show 
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an abrupt drop in lifting pressure across the tip Mach lines, the 
residual lift over most of the enclosed area being almost negli- 
gible. The attendant change in center-of— pressure location with 
Mach number may present a serious stability problem. 

As a result of this development, a wing with tips normal to the 
stream is also investigated. This case is relatively simple to 
calculate and the necessary formulas are given in closed form. 
Calculations are presented for such a wing, as well as for two 
conventional swept— back wings with 63 ° sweep of the leading edge, 
at a Mach number of 1.5. 


INTRODUCTION 

Although a number of treatments of the problem of the lift on 
a thin flat surface in supersonic flow have recently been published 
(see bibliography of reference 1, and references 2 and 3 ) the lift 
distribution over a wing with both leading and trailing edges swept 
behind their respective Mach lines has not, at the present writing, 
been determined. The only explicitly formulated solutions (all 
based on the linearized form of the flow equations) are essentially 
of the so-called "conical" type, introduced by Busemann (reference 4); 
that is, solutions in which such quantities as the velocities and 
pressure are constant along lines radiating from a single point. 

The limitation is therefore automatically imposed that the boundary 
conditions to be satisfied must also be constant along such lines. 1 
Under certain circumstances, conical fields may be superposed to 
give surfaces of nonconical plan form. Figure 1(a) shows one such 
case, a finite— span trapezoidal wing swept only slightly compared 
to the Mach lines. The basic solution here is that for a flat, 
symmetrical unyawed triangle at an angle of attack, and is conical 
with respect to the apex 0. The triangle may extend to infinity. 

The pressure distribution over the shaded portion is constant; the 
values vary between the Mach lines from 0, but are constant along 
any ray drawn through 0. In order to obtain a finite wing, it 
is necessary to cancel the pressure beyond the tips AB and A*B’ by 
superposing negatively loaded triangular surfaces with apexes at 
A and A* and one edge parallel to the stream. The conditions to be 
satisfied by the supplementary solutions are that (l) the pressure 
be constant over the subtracted surfaces in order to cancel the 
constant pressure in the basic solution, and (2) the downwash induced 
inboard of the streamwise edge be zero in order that the surface may 
remain flat. Both of these conditions are conical with respect to 
A or A’, so that a conical solution may be used. The area behind 


1 The source— distribution methods of references 2 and 3* however, are 
applicable to curvilinear plan forms within certain restrictions. 
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BOB' may, of course, Be subtracted without affecting the wing Itself 
in any way, Bince the wing lies entirely ahead of the Mach lines 
originating from any point on the surface BCB*. Consider, however, 
the case shown in figure 1(h). Here the Mach number is such that 
the Mach lines from A intersect the tips of the wing. Starting 
with the same basic triangle, we find that the area to be removed 
outboard of A includes, in this case, a region over which the 
pressure varies, and is conical with respect to 0. Since the 
boundaries of the region are, however, conical with respect to A, 
no one conical solution can satisfy the conditions. This case and 
the preceding one can be treated by the point-source method of 
reference 2. 

It has became customary to describe a leading or trailing edge 
briefly as subsonic or supersonic accordingly as the component of 
free-stream velocity normal to it is subsonic or supersonic, that is, 
as its angle of sweep is greater or less than the sweep angle of the 
Mach lines. Thus wings 1(a) and 1(b) are said to have supersonic 
leading and trailing edges. In 1(c) is shown a wing with subsonic 
leading edge and supersonic trailing edge. In this case, as in l(b), 
the pressure over the region outboard of the tip varies conically 
from the apex 0j whereas the area involved is conical from A. The 
wing shown in 1(d), with which this paper is chiefly concerned, 
presents the same difficulties, not only at the tips, but also along 
its subsonic trailing edge. 

A still higher degree of sweep relative to the Mach lines, or 
a higher aspect ratio, would result in the trailing— edge Mach lines 
intersecting the leading edge. This case, which introduces new 
problems, will not be covered in the present paper. 

The treatment of the problem will be along the usual lines of 
the linearized theories. Since the are have been so many careful 
examinations of the fundamental problem and equations of supersonic 
flow, these will be given only cursory mention where necessary to 
establish the course of the argument. The object of this paper will 
be to give a practical method for approximating the theoretical lift 
distribution for the conventional swept— back wing, and the emphasis 
will accordingly be on final formulas, where they can be given. 

The wing 1(d) with both leading and trailing edges swept behind their 
respective Mach lines, for reasons to be discussed later, can be 
treated only by successive approximations. The method of carrying 
out the successive steps will be indicated, but the amount of work 
required to find the pressure distribution beyond the first approxi- 
mation appears prohibitive without the aid of a mechanical computer. 
The total lift and moment can be found to very good accuracy 
without much difficulty, however. Mathods of approximating the 
effects beyond the first step will also be suggested. 
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The degree of approximation involved in the proposed estimates 
will depend on the aspect ratio and taper ratio of the wing, 
numerical e xam ples will he included in the paper and will serve to 
indicate the possible magnitude of the errors. 

The method of canceling the pressure at the wing tips was orig- 
inally suggested to the author by P. A. Lagerstrom of the California 
Institute of Technology. The treatment of the subsonic trailing 
edge is largely an extension of the same method of attack. Applica- 
tion of the methods to other than conventional swept wings is 
discussed in the present paper. 


COORDINATES AHD BASIC PARAMETERS 

In accordance with the usual linearized approach, the boundary 
conditions are satisfied in the horizontal plane, rather than in the 
actual plane of the wing. The coordinate system will be chosen so 
that the apex of the leading edge is at the origin, the positive 
X-axis extending back in the stream direction along the axis of 
symmetry of the wing and the 1 Y— axis perpendicular to the X— axis in 
the plane of the wing. The Z-axis is vertical, positive upward. 

If the corresponding perturbation velocities are u, v, and w 
in the same order, the pressure difference, or local lift, is then, 
in the linearized theory. 


Ap * 2pVu ( 1 ) 


4E . 4 3 * 
q. v 


where p is the density of the air, q the dynamic pressure, a nd 
Y the velocity of the stream. 

As previously noted, the conical solutions give values of the 
velocity u that are constant along. radial lines. Therefore, the 
fundamental geometric quantities are the slopes of various radial 
lines. All such slopes are measured from the stream direction as 
the reference line. It is further convenient to relate these 
quantities to the inclination l/$ of the Mach lines in each problem, 
since it is by now well known that the solutions of the linearized 
supersonic flow equation are essentially functions of Such a ratio. 
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Thus the following parameters and variables are defined: 


m 


Sigjg = 3 cot A 

slope of Mach lines 


m-fc = p x slope of trailing edge 


a = p x slope 
apex of 


of any ray from the origin, 
the wing 


or 


(3 = ^/m 2 — 1, where M is the free— stream Mach number 


Other symbols referring to relative angular quantities will be 
defined as needed in the same way. A summary of the symbols will 
be found in Appendix A. 

I - BASIC LEFT DISTRIBUTION 


WING WITH SUPERSONIC LEADING EDGE 
Pressure Distribution Over Triangular Wing 

The general procedure, as has been indicated, is to start with 
the solution for the infinite triangle having the same sweep of the 
leading edge as has the wing under investigation. This basic 
solution has only one parameter, the sweep of the leading edge 
relative to the Mach lines. Thus the incremental streamwise velocity 
distribution u^(a) over a lifting triangular plate with supersonic 
leading edges is given (reference 5) by 


MM = _-2£. cos-1 /5 5 

?a «P|/m2-I V m 2 — a 2 


( 2 ) 
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where a is the angle of attack in radians. This function is. 
plotted in figure 2. 


Total Basic Lift 

The integral of the pressure difference represented hy this 
velocity distribution over the plan form of the wing under considera- 
tion will be called the basic lift %. The integration will be 
carried out for case 1(b), in 'which the leading— edge Mach lines cut 
across the tips. The element of area is a triangle formed by two 
rays of inclination a and a + da, and haring as a third side 
either (l) the trailing edge of the wing or (2) the wing tip, 
accordingly, as a is less than or greater than at, the value 
corresponding to the ray through the tip of the trailing edge 
(fig. 2). Therefore 


Lo = 8 1" r at mt 2 co^ JJA 
qa. 3 l J 0 2 (mf-a) 2 Va a 



3 g s g ua 
2a g Ya 


da 


+ 


[* P gs2 
Ji 2a g 


Hi da 
Ya J 


(3) 


where Co is the wing chord in the plane of symmetry, ui is the 
value of va for a=l and s is the semispan of the wing. From 
the geometry of the wing we obtain 


at 


3s 

Co+ps/mt 


00 


or 

mtc D “ §§ (mt-at) 


The integration of equation (3) yields 
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Xk> _ 4ms 2 I ®t 

ZZ pSESSEE J 

qa Jt/m 2 — 1 l Lat 2 



(mfr-at) 2 ] x 1+mat 
at 2 (m-t+m) J m+at 


+ 


mt_ _ l _ (m-t-at) 2 GOS _ 1 1-mat 

at 2 m at 2 (mf-m) J im-at 


+ 2mt (mf-a-(. ) s /m2— 1 
at 2 (mt 2 ~m? Vmt 2 — 1 



mt—at 


— cos -1 



+ £ cog-lit 

m t a t 2 (mt 2 -® 2 ) rnj 


(5) 


This form becomes indeterminate when mt=aa, that is, for an 
untapered wing. For this case 

Lo = 4s g f (m+at) 2 cos -i l+™t 
qa itat 2 i^-l 4 2 m+at 


+ (m-at) 


m+at m (1-mat) 

” ^r~ 


COS - 1 


1-mat 

m-at 



+ O^t) 2 

V'mS-l 


(b L 


1-at 2 


m-at 



(6) 
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Basic Moment 

The center of pressure for this component of the lift is 
readily obtained when it is considered that each triangular element 
is uniformly loaded and therefore has its center of pressure two- 
thirds of the distance along the median hack from its apex 0. For 
a < at, the x coordinate of the center of pressure is 

2mtc 0 _ 2Ps (mt-at) 

'3 (mt-a) 3 &tfat-a) 


For a > at, the corresponding quantity is j Then the basic 

pitching moment about the apex of the wing is given by 


Mo = =3 
qa p 


L 0 0 


^JBLSfi^da + f X 3!|!^ da+ T m P!§!iiida 
3 (mf-a) 3 Ta J at 3a 3 Va J. 3a 3 Va 


( 7 ) 


or 


Mq — 8mps 3 f if at" \( _i 1+mat 

)(cos 


, l-mat 

+ cos A 

m+at mr-at 


(mt-at ) ! 


1 cos -i 1+ma t + 1 

L (mt+m ) 2 m-nat (mt-m ) 2 


cos 


— i 1— mat 

im-at J 


+ (mt-at ) 3 


f mt 2 +m? 1_ 

L (mtS-m 2 ) 2 mt 2 - 


cos - 1 i 
m 


(mt-at ) 3 v /m 2 -l f m^+m 2 mt 2 \f , 1-mtat , 1 \ 

+ ==i — - — - + — - — )\ cos x cos A — j 

(mt 2 — m 2 j/mt 2 — 1 Vm t 2 - m ^ “t 2 -! ' v mt-at mt / 


(mf-at ) g /m g — ! 
(mt 2 -l) (mt 2 — m 2 ) 


nWI 


-—at 


(mt-at) 


at 3 |/m 2 — 1 
m 2 


cosh 1 — [■ 

at J 


(8) 
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For the untapered wing 

Mo = -4mB3 3 / 3 (g^-at 2 ) (m*at) c03 _i d+mat 

<J. a 3«a-t 3 /^— 1 t m+at 


+ (m-at) 3 


5 + 3at 
4m 2 m(m-at) 2 


w ?+2 
2 (m 2 — l) 2 


cos -1 


1— mat 

m-at 


+ (m-at) 3 


m 2 +2 

-2(i^-l) 2 


_ 2 _ 

m 2 - 


cos - 1 


1 

m 


+ cosh- 1 

m 2 


1 _ _ 3 (m-at) 3 

a t 2 (m 2 — 1 ) 3 / s 


m-at / l-at 2 
2 / m 2 -! 


1 


, (2m g +l) (m-at) 

m(m?— l) 



(9) 


WING WITH SUBSONIC LEADING EDGE 

Pressure Distribution Over Triangular Wing 

The solution for the flat lifting triangle with subsonic 
leading edges has been given by Stewart (reference 6) and others. 

In the notation of this paper the velocity distribution is given by 


HA (a) 

Ya 


m g x 1 
(3E (y'l— m 2 ) /m 2 -a 2 


( 10 ) 


where E (/l— m 2 ) Is the complete elliptic integral of the second- 
kind, of modulus i/l-m 2 . It should be noted that for a given plan 
form the above expression and all those that follow in th is section 
vary with Mach number only as the reciprocal of E (/l— m 2 ) . 
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Total Basic Lift 

Integration of equation (10) over the wing plan form gives 


Lo _ 4s 2 

“^(“f-at) 2 / “t 


1 a at^E (v/l-^2) 


cos — COS x 

L m (mt-at ) mt _ 




( 11 ) 


As Before , a special form is required when mt=m : 


L o _ 4 s£ f (m-at) 2 

< l a ( \Jl-nP- ) \ 3 


2m-a-t 

L (ot-at ) 2 


\/m^-at 2 — 2 


+ at 


V^-at 2 j- 
( 12 ) 


Basic Moment 


The moment is obtained from equation (j) with, however, the 
third integral deleted, since m < 1, and with m substituted 
for 1 as the upper limit of the second integral. Setting ^ 
equal to the expression given by equation (10) we obtain 


Mq _ - Um 2 ps 3 1 

(mt-at) 3 

j 2mt 2 +m 2 

COB -i me-^ntat 

qa 3 at^ ( v/l-® 2 ) j 

(m t 2-mB ) 

l(m t 2_m2) 3/2 

m (mt-at) 


— cos' 




3m 


m,. J mt 2 -® 2 mt ; 


4,2 J 


mt 


mt 2 -m? 


ipg-mtat + 3 (mt-at ) g 
m 2 mt 2 -n? 


Vm 2 -^ 


2 + coshr-i JL 
a a t 


( 13 ) 
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For the untapered wing the corresponding formula is 


Mo 8fia 3 

3at^ (/b-®*) *-3 


ij 


(»«t) + 2i + g fe ag 

10 5m 




■-at' 


_ Ifestl! + s£. OOBh -i JL 1 (14) 

15m 2m at J 

II - TREATMENT OF THE TIP 

In canceling the pressures at and "beyond the tip, the problem 
is brought within the limitations of the conical solutions by 
breaking the region down into an infinite number of constantly loaded 
overlapping sectors of infinite extent. (See fig. 2.) These sectors 
are bounded on one Bide by the wing tip; the second side is the 
extension of a ray from the apex 0 of the wing. The (constant). 

pressure on each sector is then 2pY— ^da, where a is the 

da 

previously defined parameter 0 times the inclination of the ray 
from 0, the ray in this case being the particular one bounding 
the element. 

Any element bounded by a ray of inclination a/0 will hereinafter 
be referred to as the element a, or sector a. The apex of the 
sector will be designated by the coordinates x^, y a and the value of 
the constant incremental velocity over the element by 1 %. 


ELEMENTARY SOLUTION FOE TIP 

The flow fields to be superposed on the basic lifting triangle 
to cancel lift beyond the tips will be required to satisfy the 
following conditions: 

1. The streamvd.se velocity u must be constant over the 
sector between the tip and the specified ray. 

2. The velocity u must be equal in value but opposite in 
sign on the upper and lower surfaces of the sector, in order to 
produce lift. 

3. The streamwise velocity u must be zero outside the sector 
on the free stream side. 
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4. The associated downwash or vertical velocity w must be 
zero outside the sector on that side adjoining the wing, in order 
not to disturb the boundary condition of constant slope of the wing 
already satisfied by the basic solution. 

5. All induced velocities must go to zero on the Mach cone 
from the apex of the sector. For a discussion of this condition, 
see reference 7. 

6. The velocity components .u, v, and w must satisfy the 
linearized supersonic flow equation and the equations for irrota- 
tionality. 

The solution of the supersonic flow equations satisfying the 
specified boundary conditions (l through 5) has been derived by 
Lagerstrom, using the method of analytic extension. No attempt will 
be made to reproduce the derivation here. The discussion which 
follows is included only to illustrate the nature of the solution. 

Condition (6) may be shown to be satisfied by any analytic 
function of the complex variable 


y+iz 

x+Zx 2 — P 2 y 2 -P 2 z 2 


(15) 


It is convenient to set up the solution for u as such a function 
of e, since the boundary conditions are largely specified in 
terms of u. Condition k, a condition on w, can also be related 
to u by the following procedure: 

Let u represent the real part of the complex analytic function 
U(e) = u(x,y,z) + iu(x,y,z) and let the corresponding complex poten- 
tial be 


$ = J' U dx = (p+icp 

Then w = Bcp/dz will be the real part of W = d$/dz. We may write 


as 


*-fc/ 

1? 


U dx 


dx 


(16) 
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the additio nal terms arising in the differentiation under the 
integral sign "being zero. 

The lower limit of the integral is arbitrary as far as the 
determination of $ is concerned, since it can only change $• by 
a constant and its derivatives not at all. However, in order that 
condition 5 be satisfied it now appears necessary to specify a 
point on the Ifech cone as the limit. 

Equation (l6) may now be written 


de 

Be fbx. 


The ratio of partial derivatives appearing under the integral 
sign is independent of the function assumed for U. In the plane 
z=0, it is equal to 



(See footnote 2 . ) 


Q£ / (Jn 

Inside the Mach cone x > |3y and is an imaginary quantity. 

de /ox 

It follows that over any region of the X— Y plane In which U has 
no imaginary part, the real part of W will remain constant. Since 
the region to which condition ^ applies is bounded by the Mach cone, 
on which W is zero, it follows that condition H will be satisfied 
if the imaginary part of U Is zero over the region. 

It is convenient now to introduce the variable 



which will be seen to reduce to a =■ 0 ? when z=0. 


2 The general expression for 


Be/Bz 

Be/Bx 


indeterminate, rather than zero. 


includes a real part which becomes 
when both y and z equal zero. 


and the argument that follows no longer holds. However, the X— axis. 


along which this occurs, is not included in the region under dis- 
cussion (see condition h), so that the reasoning is adequate. 
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Without rewriting equation (15 ), we will consider the variable 
e to he referred to the apex (x a ,y a ) °f the elementary sector a 
as origin and define the new angular parameter (corresponding to a 
in the fixed reference coordinates) 


t a = 3 fp 1 ( 18 ) 

x— x a 

This will then he the value of £ in the plane of the wing. The 

ray t a associated with any point P(x,y) and a particular sector 

a is shown in figure 2. When the solution is being applied to the 

region of the wing tip, y a = s and x a = — S o that 

8 . 




x—ps/a 


(19) 


In the region of the right-hand tip, t a is negative on the 
wing, and positive on the elementary sector. It takes on the 
following special values: 

t a = a on the free— stream boundary of the sector 

t a = o along the wing tip 

t a = +1 along either Mach line 

Conditions 1 through 4 may now be written 


1 

2 



3. 

4. 


o < t a < a, lim u(x,y,z) = — 1% (1% constant) 
z-> 0 | z l 

a < t a < 1, u (x,y,o) = 0 

,-l < t a < o, u (x,y,o) = 0 

Consider novr the complex function 


U = ^ cos" 1 a+ ^ +2a ^ = % 1 1 a+j^+2a^+2|/a^ (l+a) (l+£) 

* ^ * 8 — 


( 20 ) 


As z approaches zero through positive and negative values, £ 

can be seen to approach the real value t a through conjugate complex 
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values. The argument of the logarithm in equation (20) will also 
take on conjugate complex values, as will the logarithm itself. The 
real part of TJ, however, is proportional to the imaginary part of 
the logarithm and will therefore differ in sign above and below the 
X-Y plane, as required by condition 2. In the plane z=0. 


u = ± r.p. 


u a 

T 


cos -1 


a+t a + 2 at a 


(< 51 ) 


The arg ume nt of the inverse cosine in equation. (21) is plotted 
against t a in figure 3 . When the argument is greater than 1 or 
less than - 1 the function JJ is complex. Both the real part 
u and the imaginary part u are plotted, but only for the upper 
sign in equation ( 2 l). 

It will be observed that conditions 1, 2, and 3 are satisfied by 
u and condition k by the corresponding imaginary function u. By 
equation (l 6 ) and the corresponding expression for the lateral 
velocity component, w and v are determined, regardless of the 
function assumed for u, in such a way as to satisfy condition 5» 

That u also satisfies that c ondit ion is readily demonstrated, 
since on the Mach cone x = Pl/y 2 +z 2 and £ - x/Py > 1. When £ is 

greater than 1 , the argument £ ie positive and greater than 

£■— a 

one, so that U (equation ( 20 )) is imaginary, or u= 0 . 

APPLICATION TO WING WITH SUPERSONIC LEADING EDGE 

As long as the wing tapers in the usual sense, that is, from 
root to tip, the condition of a supersonic leading edge implies a 
supersonic trailing edge as well. In that case the pressure in the 
wake may be canceled without inducing any pressures on the wing ahead. 
The section that follows will therefore cover completely all correc- 
tions to the basic lift for a swep1>-back wing of type 1(b). Moreover 
when the total lift has been found, the drag due to lift is also 
known in this case, since it is merely the lift times the angle of 
attack when the leading edge of the wing is supersonic. 


Induced Pressure Distribution 

The basic lift distribution has been given by equation (2). 
Following the procedure outlined in the preceding section, we divide 
the lift outboard of the wing tip along projections of rays from the 
apex of the wing (fig. 2 ). Beyond a=l no further division is 
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necessary, since the basic pressure is constant at ux ahead of the 
Mach line from 0. Thus the solution to be superposed may be con- 
sidered to be made up of the sector a=m with constant incremental 
velocity Ua = — u x and an infinite number of narrower sectors 
(a<l) each with constant infinitesimal velocity 


du& 
= Ta 


da 


( 22 ) 


Then the total induced correction to the basic velocity distribution 
is S 




where t a is the function of x,y and a defined by equation (19), 
and must be so expressed before integrating; and where tm is the 
function of x and y obtained by setting a equal to m in t a . 
The lower limit a© of the integral corresponds to the most rearward 
sector that can influence the pressure of the point x,y; that is, 
the sector such that the Mach line from its apex x a ,s passes 
through the point x,y. The' expression for a©, found by setting 
t a equal to —1 and solving for a, is 


a 0 


3s 

x+0 (y— s) 


(sit) 


This par ame ter will be additionally useful as the value of a at 
which the function given by equation (21) goes to zero and its 
derivative has a singularity. If a 0 is equal to or greater than 1, 
the induced velocity Au reduces to the first term of equation (23). 

From equation (2) 


and 


(25) 


duA _ 2mVq a 

da n0 (rn^-a 2 )\/l— a 2 



(26) 


The integration of equation (23) for the local pressure has 
not been carried out, since the solution for this case is obtained 
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mor e easily by the method of reference 2. However, the total lift 
induced on the wing by canceling pressure outboard of the tip is 
readily obtained by the present method if the order of integration 
is reversed and the total lift induced by a single element found 
first. 


- Total Induced Lift 

The change in lift due to any tip element a, obtained from 
the velocity distribution (equation (21)), is 

<OL)a = 2pV f° Sk cos-i t+Vgat a dtg ( 27 ) 

J—x t a — a dt a 


where dS is the element of wing area. If the line t a = —1 does 
not cut the opposite tip. 


dS _ frnt g s a f l l\ 2 1 

dt a 2 \at a / (mt-t a ) 2 


(28) 


where at is the value of a defined by equation (4). Then 
equation (27) becomes, after integration by parts, 

(AL) a = % (at-a) 2 g (a) (29) 

where 

Va+a g a 

v/mt+mt 2 m t 


g(a) = — i— 
m t -a 


and the total change in lift is 


(taking both wing tips into account) 


AL 



pYmt 2 3s 2 
“a * 5 



duA (at-a ) 2 
da a 2 



(30) 


For the first term of equation (30), a is taken equal to m in 
equation (29 ) , and u a equal to m (equation (25)). The second 
term, is integrable by elementary methods. The final result is 
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AL = 2m b 2s2m ^t) 2 f 1_ cos _ 1 

qa * a + 2 yKI \Mn-t-a) W 


1 ~ ma t _ 1 / m+m 2 cos _ 1 m - gma t 4a t 

m-at m V mt+mt 2 m-at 


(m+af.) 2 f 1_ coa _ x 1+mat _ 1 /jEsL 

i (mt+m) ^mt m+at m y m-^+rnt 5 


L coa -i -^a+2mat+at 
; 2 m+at 


2 ^-a^/m 2 -! 
m-t (mt 2 -® 2 Vmt 2 — 1 


m+at m y mt+mt 2 m+at 

' coa -i _ coa -! sjw s 

v mt— at mt— at 


As before, a special form is needed for the untapered wing: 


Ut> 

l a Jtat 2 /™^— 1 


(3m 2 -!) (m-a t ) 
2m (m 2 —1 ) " 


l-ma+ 
cos-i - 

m-at 


- (mnat ) 


(3mr-2) (m-at) 
2m(m-l) 


2 cos -1 


m-2mat+at 


(m-tat) 2 / cog-i 1+ma t _ F&-1 cog -i ~ m+gma t+ a t 
2m \ m+at v m+1 m+at 


(m-at)y/l— at ( /a t \A+ a t 


/m^-1 


/m+1 


Induced Moment 

The moment due to this correction to the lift may also be 
expressed in closed form, but the formulas are so long and unwieldy 
that their integration in general terms was not considered worthwhile. 
The moment arm of any element dt a from x^_,s about the Y-axis is 



The moment induced by any element a is 
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W, - u* [ «M ♦ h«] (33) 

where 

. / \ /, mt \ / \ 1 / a+a 2 

(a) - ^ + mt-ay 8 ' a ^ 2 (l+mt) (mt-aT V mt+mt 2 


The total induced moment (including both wing tips) is 


AM = -2 


+ ^ [ 8(a) + If h(a) ] & 


As in the case of the lift, the first term is evaluated hy setting 
a equal to m and Ua equal to ui in equation (33). In the case 
that mt=m, special limiting forms of g and h will he required, 
namely. 


eM 


and 

h (m,m) = — i — \l 3 — _ 

m (1+m) |_ 8 (l+m) _ 

Integration of the second term can he done numerically for any 
particular case. 


APPLICATION TO WING WITH SUBSONIC LEADING EDGE 

The case of a wing with subsonic leading edge (m<l) requires 
slightly different treatment because the pressure at the leading edge 
is infinite. The formulas to he obtained for this case will, how- 
ever, give the complete theoretical solution for such a wing as long 
as the trailing edge is not also subsonic. Where this condition does 
not hold (m^cl), the following sections will give the first step in 
the complete solution. 


Induced Pressure Distribution 

It is necessary to write equation (23) for the induced velocity 
at x,y as 
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r-UA (X) , X+t-» +2\t-], 1 P ^ du A , a+t_+2at„ I 

Au » 11m JHl_i cos" 1 A — + - / — ^ cos" 1 &*■ & da 

X ->m w *fc\— X Ja o da ‘to,— a -* 

( 35 ) 

where, as before, t a must he replaced by a function of x, y, and 
a, (equation (I 9 )). Integration of the second term by parts gives 
rise to a term which exactly cancels the first term and is zero at 
the other limit, leaving after substitution for u^. 


= m 2 (x+Bv ) feo ( B -y) , r m 1 % ( 36 ) 

Va itE (/l-m 2 ) * s Ja. Q (Py— ax)/ (m^-a 2 ) (l+a ) (a-a 0 ) 

This integral is finite and can be evaluated in terms of 
elliptic integrals as follows: 


Au 


m 


3/2 


Va JtE (/OT)p 


/ 


gP-ferZ). K(k) 

x+py 


— 2x 


/. „ P- -f F (tjk* ) £ _ K(k)E (k* ) 

Vm2 X 2-p2y2 \ K (k») 2 ' ' 1 1 


+ K(k)E(t 


,k')} 


(37) 


where the modulus k = , k* = /l— k 2 


/a 0 (mx+Py) 

and the argument i)r = sin -1 . / -AL. — , The elliptic integrals 

V Ps(a 0 +m) 

F, K, and E are extensively tabulated in reference 8. They are also 
available in Pierce's Table of Integrals and most handbooks, but 
experience has shown that k* is generally very nearly unity and 
that K(k') as tabulated in the handbooks is not always satisfactory 
for interpolation. The following series (reference 9 ) converges 
very rapidly and is preferable for computing K(k' ) when K is small: 


K(k») = -log | - i(log | -hDk 2 - ^(log ^ + £)k 4 -... 


At the tip, y=s and the first term in equation (37) vanishes. 
In the second tern, t becomes £ and E (i|r,k* ) and F (t,k' ) 

reduce to the complete integrals E (k* ) and K(k* ), respectively. 



BACA TR No. 1555 


21 


The induced Telocity correction Au will then "be seen to be exactly 
equal to — uAj bringing the lift to zero at the wing tip. 


Drop in Lift Across Tip Mach Line 


An interesting effect shows itself at the other limit of the 
tip region — the Mach line from the tip of the leading edge. Along 
this line only the influence of the leading-edge pressure is felt; 
that is, a 0 =EU Then k=0, k’=l, l(k) = — . K(k* ) = » and E (k* )=1. 


I(\|r,lc*) reduces to sin t 


/ bx+$y 
V 2gs 


and finally 


Au _ -m. s 

Ta 1 (y/l-m^ ) \J2 (x+gy) (rax-gs) 


(38) 


This result corresponds to a finite drop in pressure across the- Mach 
line from the tip, an effect which does not appear as long as the 
leading edge is ahead of the Mach lines. The magnitude of the drop 
relative to the pressure just ahead of the Mach line should be of 
interest. This ratio is 


Au _ — gs ! 6n?-a g ) 

UA m V 2 (x+gy) (mx— gs) 

/ A pS 

Since a=Py/x, and on the Mach line from the tip P (s— j) — jp 

the ratio Au/ua can be rewritten very simply in terms of m and 
a as follows : 

Au = _ / (l+a) (m+a)' /oq) 

UA V 2m (l+m) 


This function is plotted in figure 4 against a/m. 

It will be seen that the percentage of loss of lift at the tip 
is very, large and, in fact, that for any but the lowest aspect ratio 
wings the lift remaining in that region, which drops to zero at the 
tip itself, is almost negligible. Some indication of this effect ig 
contained in the results of reference 10 for the limiting case of 
m=0. It is an effect of considerable practical interest, suggesting 
as it does a considerable change in pitching moment with Mach number 
as well as the inefficiency of the tip area of the wing. 
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Total Induced Lift 

Proceeding to the calculation of the total correction to the 
lift, we integrate the change in lift (AL) a (equation (29) ) 
induced by each element a, over the range at < a < m. Again 
we encounter the difficulty that u^ -> <*> as a -£• m and must 
therefore write the total lift in terms of limiting values,- For 
convenience the function 


G(a) 



(40) 


from equation (29), is defined. Then the total induced lift, corre- 
sponding to equation (30) for the preceding case, may he written 


AL = 2pVm-t 2 (3s 2 lim 
a->m 


-u A (a)G(a) + 


/at 


G (a) da 


(41) 


Integrating hy parts results in cancelation of the first term 
inside the brackets. Since G(at) is zero, equation (4l) reduces to 

AL = -2pVmt, 2 ps 2 / u A (a)G' (a) da (42) 

J at 

Equation (42) has been integrated in terms of elliptic integrals, 
hut because the result involves several new functions it was thought 
better to present it in an appendix (Appendix B) . For practical use 
graphical or numerical integration may be preferred. Here again, 
however, the difficulty arises.. at the leading edge that u,^ and 
therefore the integrand of (42) becomes infinite, although the inte- 
gral is finite. For numerical integration it is possible to dispose 
of the singularity by integrating once more by parts. A somewhat 
simpler method is to write 

AL = — 2pVmt 2 ps 2 I f u A (a) 

5 '•■ '■■■I : ^ a t 

The integrand of the first term can now be shown to approach 
zero as a approaches m, and the second term is readily inte— 
grable, since G* (m) is a constant. The result, when u A has 
been replaced by its value from equation (10) and the second term 
integrated, is 


G* (a)-G* (m) 


da + 


r> m 

/ u A (a) G’ (m)da 
J at 
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AL -^ m ^ m t £a2 
* a _ E(v^) 


G’ (m) cos 1 — 
m 



£1 (a)-fi* ( mj da 
i/mS-a 2 


(43) 


G’(a)-2=|% [("21 + 2fc2t') g(a)_ 

»t a 2 L\a mt-ay 


±±t_ 


2 (mt-a)/ (mt+mt 2 ) (a+a 2 ) J 


m 


When n%=m, G* (m) reduces to 


m-at 

a^m 4 (1+m) 


1 + 5±t 

8a^ (1 -hd.)_ 


Induced Moment 

The moment Induced hy cancellation of the tip pressure is found, 
as in the preceding case, hy the integration of (AM) a (equation (33)) 
oyer the appropriate range of a. Again, the actual Integration is 
test done numerically for any specific case in which the moment may 
he required. Because u a becomes infinite at the leading edge. It 
will first he necessary to Integrate hy parts, as in the case of 
■the total lift. Following the same procedure, we let 



g(a) 


a ~at 

3at 


h (a) 


(45) 


from equation (33 ) > 

AM = 2pVP®s ^mt^lim 
a->m 

= — 2pV0 2 s 3 n% 2 


so that the total induced moment is 


' a du^ 


-u A (a) H(a) + P H(a)da 

* Ja t da 



(46) 


For numerical integration, the procedure outlined for finding the 
lift is repeated, giving 


AM -4 m?m t as3 P 
i a E (\/i^) 


H* (m) cos 1 ~ 



H* (a)— H* (m) 
y/m^-a 2 



(47) 
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H' (a) 


a-at 

a t e a 3 


2at a— at , gt (a— at) 2 

a m t -a 3&t M) 2 


6 (a) 


4 2=24 (2at 4 £fc£i) h (a) 
3 \ a mt-a/ 


a-a t 


2 (mt-a)\/mt4-mt 2 \/a-*-a £ 


l 3a t 


1 4 - 


m t 


mt-a 2 (l+mt) 


When mt=m, H* (m) becomes 


ifr-at f at 

at 2 m 4 (1+m) { m 

V 


(m-at ) (X +2m) 
6at (14-m) 


1 4 at 

2 m 


(m-at) (l+2 m) 
6m (l+m) 



2nr+ -i2g!L~ 
m (l+m.) „ 



III - WING WITH SUBSONIC TRAILING EDGE 

The problem of satisfying the Kutta condition of zero pressure 
difference at the trailing edge and in the wake behind a highly 
swept wing is attacked in a manner similar to that used, on the tip. 
The region behind the actual wing surface is divided into segments 
over which constant pressure increments may be assumed by extending 
the rays from the apex of the wing through points Xa,y a on the 
trailing edge. This division of the wake region is indicated in 
figure 5j with a typical elementary sector a indicated by heavy 
lines. 

ELEMENTARY SOLUTIONS FOR TRAILING EDGE 

The element in this case differs from that in the preceding 
one in that the fixed side is no longer parallel to the stream but 
is inclined at the angle tan -1 ^ to it. The boundary conditions 

to be satisfied by the - elementary solution are therefore (on the 
right-hand half of the wing) 
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} 


3. 

k. 


cd 

-P 

VI 

cd 

< 

“t 

u 

= ±% 1 


— 1 < t a 

< 

a 

u 

= 0 

1 z — Of t - t a 

“t <t a 

< 

+ 1 

u 

= 0 



in addition to the general conditions 5 and 6 stated previously. 
The pertinent results of the preceding discussion will he assumed 
in order to shorten the discussion at this point. 


The variable t a is that defined by equation (l8) and refers 
to the slope of rays originating at the apexes of the trai ling-edge 
elements a. Along the trailing edge 


y a 


2 k 

P 


fra-Co) 


Since 


D ya 

a=P- — , we 
X a 

xaay 

solve 

for 

: constants 


and 

Co: 




x a 


_ m b c o 
mf-a 




Then 


Py a = 

“t-a 


gy (mt-a) - mtc Q a 

x(mt-a) - mt c o 


as functions of a 


m 




The solution for the oblique triangle satisfying the specified 
conditions is (in the plane of the wing) 


u = r.p. cos-i ~ feglik^a) (50) 

77 (l-mj.) (t a — a) 


Equation (21) will be seen to be the special case of equation (50) 
in which m-^0, with the axes reversed because the relative positions 
of wing and element are now reversed. (See fig. 5») 
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At this point a difficulty not encountered in the treatment of 
the tip must he considered. This difficulty is most apparent in 
connection with the element a=0, that is, the right triangle with 
apex at the trailing edge of the root section. The Mach cone from 
that point includes segments of both wings, so that condition 4 should 
properly be applied over part of the negative range of t a as well 
as the positive. This could be done for a=0, but for any other 
small value of a, the area to be so treated would be nonconical 
with respect to the apex of the element a and therefore the bound- 
ary condition of zero downwash could not be fitted into the conical 
solution. Thus a deviation from the flat— plate boundary condition 
seems unavoidable. However, the error can be minimized by the choice 
of a solution for a=0 that does not introduce any downwash on 
the far wing. 

The pressure at the root of the trailing edge is most simply 
disposed of by means of a symmetrical solution satisfying the 
following conditions: 


-mt <t Q < + mt 


u — u 0 ■ (o) 


-1 < to < -mt 

+mt < t D ^ + 1 


u = o 


where t 0 refers to rays emanating from the trailing edge of the 
root section (a=0, Xa=co, ya=0) . The required solution is 


u 


r.p. Ho cos-i .(l + Pt 2 )tof-2m t g 

(l-m t 2 )t 0 2 


(51) 


From equation (10) 


u 0 = — (52) 

0E (/l-m?) 

The major portion of the pressure in the wake will be canceled 
by solution (51) ; the remaining pressure may be canceled by the 
solutions given in equation (50) where u a is, as before, da. 

The downwash induced on the far wing by the solutions should be, and 
has in fact been found to be, small, first, because du^/da is small; 
and second, because as a and 1% increase , the portion of the far 
wing surface affected by the element a decreases in area and 
increases in distance from the element. 
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APPLICATION AT TRAILING EDGE 
Pressure Distribution Hear Trailing EdgO 

In determining the pressure at any point x,y on the wing 
surface it is first necessary to determine the most rearward 
element a Q that will influence the point x,y. As before the 
value of a D is found by setting t a (equation (4 9) ) equal, in 
this case, to 1, so that 


a Q = mt Q — (53) 

Py-Wc^t “ x 


The velocity induced at any point within the Mach cone from the 
vertex of the trailing edge by superposition of the trailing— edge 
solutions discussed thus far is given by equation (51) plus the 
integral of equation (50) with respect to a, from a=0 to a=a D . 
The total is to be subtracted from the basic distribution given by 
equation (10) . 


Total Induced Lift 

If the Mach lines from the apex of the element a do not cross 
the leading edge of the wing, the conical element of area for the 
integration of the pressure to the tip is 


dS = P m t 2 s g (at -a) 2 dt a 
2a+ 2 (mt-a) 2 t a 2 


( 5 ^) 


For the symmetrical solution, a=0 and 


dS _ 0s 2 
dt 0 2t 0 e 


Integrating equation (51) over dS from t 0 =mt to t D =l and 
multiplying by two to take account of the negative half of the 
wing gives, for the symmetrical solution. 


£L = — 4s 2 m 

<l a mtE (/l-m. 2 ) 



( 55 ) 
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When nit ia nearly liO, the expression In the brackets may he 
difficult to evaluate accurately. The following series expansion 
(from reference 8) is useful: 


1 (/1-mt 2 ) = 

JC c. 




+ 


i g 3 2 

224252 


E 4 


+. 



where 


E 

1 +mt 


The total correction to the lift induced by one of the oblique 
triangles in the wake is found by integrating equation (50) over dS 
as given by equation (54), from the trailing edge to t a =l, The 
resulting formula is 




— pVmt g s g p (at-a) g ua 

at 2 (mt— a) 2 




(mt-a) (l— a) 
mt 



For the total induced correction to the lift due to the oblique 
elements, (AL) a will have to be integrated over both wings, after 
substitution for 1%, and added to AL obtained from equation (55)* 


The integration has been done for both tapered and untapered 
wings and results in the following formulas: 


For mj. / m, 

AL — 4m t s g 

qa at 2 (mt 2 — m 2 )E (y/l— m 2 ) 


fflg(m t -a t ) 2 

\J mt 2 — in? 


coe-i - coa-v JL 

m (mt— at ) - 


+ 


+ 


(m 2 -m t a t) • |2m-ym2-a t 2 1 + | - 

l/“t(l-m) K+m) {i ( m+at ) g 2t^ + (m-^) 2 liS 

L m^+m 1— m . 


K-at) 

F(® t , k) 


F ($b,k) 
1 


“t 




(m 2 -m t a t ) £ + 3m 2 (114.-^ ) s 


E($ t ,k) -E($ 0 ,k) 


} 

(57) 
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/mt +m „ , / 2m(l-mt) 

,/— — and k = / -7 r-r r 

V 2]% Y ( m t +m ) 


m 2 ~a t 2 ^ (m-a^) 2 ^ 5 __i_ \ 

~2n^ + ~ 8m 2 Y° l- m / 



(ni-a t ) 2 (5m-3) ' 
16m 2 (1-mT - 


3-m 

1+m 


(58) 
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Induced Moment 


Any constantly loaded element of area 
has a moment arm about the Y— axis given by 



from Xa,y a 


Bs 

a t K-a) 


r K-t) * 


( 59 ) 


For the symmetrical solution, a =0 and this expression reduces to 


Ps 


4 _ 2 _\ 

a t“t 3t cJ 


The moment induced on both halves of the wing by superposition of 
the symmetrical solution is, from equation (51) 


& , 2pTfis ’' 0 r 1 

/ mt-at 

* Jmt 

V. at^t 

and 


AM ^ 4 mgs 3 

mt-at 

qa (\/l-m 2 ) 

a t 



cos— 1 


(l+mt 2 ) to 2 — 2mt 2 
(l-mt 8 )t 0 2 


1 - |eVi^) 


+ 


1-mt 2 

~T~ 



(60) 


The mon»nt induced by a single oblique element a is most 
conveniently written in two parts, corresponding to the two terms 
of equation (59): 


(£M) 


a,i 


-gs (m-fc-afr) 
at (mt-a) 


H 


(61(a)) 
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In (AM) a ^ 2 ^ the corresponding bracketed expression equals 


1-mt 

lint 2 



1+mt mt\ 
~ + a ) 



Q + 


1.3.5 

6.8 



SUCCESSIVE STEPS FOE COMPLETE SOLUTION 

Application of the procedure described so far will result in 
cancellation of the pressure along the forward part of the tip and 
along the inboard portion of the trailing edge. Over that part of 
the tip that lies within the Mach cone from the apex of the trailing 
edge, however, the method of removal of the trailing-edge pressures 
has induced an extraneous negative pressure, for which an additional 
correction will have to be made. Similarly, the tip solutions have 
introduced negative pressure along that portion of the trailing edge 
falling within the Mach cones from the tip. There is also a complex 
distribution of pressure in the wake and in the stream beyond the 
wing hips. 

The lift restored to the wing as a result of the cancellation 
of these erroneous negative pressures will be referred to as secondary 
corrections to the basic lift. 

Figure 6 shows the regions near the wing tip affected by these 
secondary corrections. The errors introduced by the primary correc- 
tions occur along and behind the trailing edge from t m = —1 to the 
tip, and along and outboard of the tip from t 0 = 1 to the trailing 
edge. Since in general very. little of the wing surface lies within 
the Mach cones (shown shaded in fig. 6) from these regions, and the 
actual flow pattern of the viscous fluid would differ considerably 
from what we would calculate, it is questionable hew far one should 
go in attempting to correct for the residual pressures there. The 
exact procedure will be outlined here, and may be feasible when a 
mechanical computer is available to perform the integrations 
numerically. A method for approximating the pressure distribu- 
tion, lift and moment will be suggested in a following section. 

It should first be noted that the distribution of residual 
lift in the field is no longer conical with respect to any one 
point, but is composed of an infinite number of superimposed conical 
fields originating at various points along the trailing edge and tip. 
In order to cancel these pressures accurately, it is necessary to go 
back to the application of the elementary solutions at the wing tip 
and trailing edge and cancel the extraneous pressures introduced by 
each one. This will involve, for each element, the integration of 
an infinite number of elementary solutions of the other type; that 
is, trailing-edge— induced pressures will be canceled by the applica- 
tion of tip solutions, and vice versa. 
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Consider an element at the trailing edge (fig. 7) with uniform 
pressure Uq,, giving rise to a pressure field defined hy equation (50) 
in terms of t a , the angular displacement around Xg^ya on the 
trailing edge. The velocity 


(to). - ^ ooa-l (ta-Qt)-(mt-o) O-ta) 

a * (1-^) (t a -a) 


(63) 


with 


n y fc _y a 


(64) 


induced at points x^y^ on the wing tip hy the individual element 
a may then he removed hy superposing triangular elements at the 
tip formed hy extensions of the rays t a from Xa,y a on the 
trailing edge. Such tip elements in turn induce velocities within 
their Mach cones which may he expressed in terms of the angular 
displacement around the apexes of the element *h , y v This dis- 
placement may he denoted hy the parameter 


% 


= P 


y-yt 

X-Xh 


(65) 


which is a function of a and t a * Then, from equation (21), the 
effect of a single tip elenent hounded hy the ray t a is proportional 
to 


1 cos -i V^-ath (66) 

* t b -t a 


Since ^u) a =0 at t a = 1, the total velocity induced at a point 
x,y on the wing hy removal of all tip pressures resulting from the 
wake element a will he 


p ^h 1 d (Au) a t a +tb+2t a t"b 

— 1 / cos t — t dt a 

* v/t R=1 dt a tb-^a 


(67) 


The derivative d(Au) a /dt a is infinite at the lower limit. 

In order to perform the indicated integration graphically it is 

preferable to write - ^ ^ dt a * d(^u) a and Integrate hy 

at a 
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plotting the inverse cosine function against (Au) a for the 
indicated range of values of t a . 

The foregoing procedure must in turn he followed for values of 
a from a=0 back to that value of a for which t-^ (t a ,x,y) = —1 
when t a =l (i.e., for which the point x,y lies on the reflected 
Mach line from x^ya)* and the results integrated with respect 
to a. 


The integrated expressions, ( 29 ) for the lift and (33) for 
the moment due to each tip element, apply to the elements from 
xjj,yb if a is replaced in the functions g and h by t a and 
elsewhere by a (x^y^) = Py^/x^. In the coefficient, — d (^u) a /dt a 
will replace Ua as the constant velocity over the section. For 
the total increment in lift and moment caused by correcting for a 
single wake element a, and (^M)h obtained in this way 

would have to be integrated with respect to t a from t a = n% to 1. 
Again, since d (^u) a /dt a is infinite at the limits, it is prefer- 
able to integrate with respect to (Au) a . The results would have 
to be integrated again with respect to a, from a=0 to a=a-(-. 

In canceling the tip pressures introduced by the symmetrical 
solution (equation (51) ), the same procedure would be used with 
(Au) a replaced by 


uo , 

COB * 

Jt 


(l+mt 2 )t 0 e -2m t 2 

(l-m-t^Ho 2 


and t a = t D , x a = c 0 and y a = 0. Integration with respect to 
a is, of course, not required. 

A similar procedure to that Just outlined would be followed 
in canceling the extraneous pressures introduced at the trailing 
edge by the original application of canceling solutions at the 
tip. In this case, however, the presence in the basic solution 
of infinite velocities at the leading edge leads to considerable 
difficulty in obtaining a satisfactory method for numerical 
integrations. Since the results presented in figure !)■ suggest 
that the wing tip should be raked inward for efficient design, 
it does not appear worthwhile to develop this part of the solution 
in any further detail. The approximate method for treating the 
extraneous pressures at the tip presented in the following section 
is considered satisfactory for practical investigation of the 
conventional wing tip. 

Where the Mach number is very low (close to 1.0), it might be 
necessary to carry the foregoing procedure through another step. 
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canceling the pressures introduced at the opposite edge by the new 
solutions in t^. It is clear that the process outlined is a 
converging one, since (l) the functions giving the pressure distri- 
bution are everywhere smaller in value than the pressure being 
removed, and (2) the area affected by each successive step is 
smaller than the preceding. The approximations suggested in the 
following section may be made at any later step in the calculations. 


METHODS FOB APPROXIMATE SOLUTION 

The following method of approximating the secondary corrections 
will result in cancellation of the extraneous pressures at the 
trailing edge and along the tip of the wing, but will not take 
account of the exact variation of pressure in the stream. It will 
be apparent however, that the residual errors will be small, partic- 
ularly in their effect on the wing. 


Correction for Tip-Induced Pressures 


The assumption will be made that since the pressures induced 
by the tip are in the main due to cancellation of the infinite 
pressure at the leading edge, the extraneous pressure field introduced 
by the tip solutions will be very nearly canceled by a pressure field 
conical from the tip of the leading edge (fig. 8), provided the 
pressures are made to cancel exactly along the trailing edge. In the 
conical pressure field to be superimposed, the velocities are held 
constant at Au (equation (37) ) along the rays tm projected from 
the tip of the leading edge back into the wake. To effect this 
approximate cancellation, values of Au are calculated from 
equation (37) for points xb,yb along the trailing edge. The 
corresponding values of % are 


tja 


yb - 0 


x b 


_ M 


m 


(68) 


1 / 3 1 the particular point at which the line tjn = —1 intersects 
the trailing edge be designated by x*, y* and other symbols 
referring to that point be similarly starred. Then the velocity 
induced at any point x,y on the wing by removal of Au along the 
trailing edge will be (from equation (50) ) 
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-Au* 2(t*-m t )-(m t +l) (l_t*) 

~ cos 
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dAu COB — 1 ^ ^m) )~ (l— ^b) 

dt m (V^m) 


(69) 


where Au* is given hy equation (38), t* is the parameter 


and t^ is the previously defined parameter (equation (65)) measuring 
displacement around other points x^y^, on the trailing edge. The 
upper limit of integration is a function of the point x,y at which 
the pressure is being determined. 

The derivative dAu /dtm would have to be determined numerically 
or graphically from a plot of the calculated values of Au against 
tm» In order to avoid this procedure, it is preferable to rewrite 
equation (69) as 


-Au* _ x 2(t*-^)-K+l) (1-t*) 
* 008 (1-ffl-t) (t*+l) 


_ i f 1 cos-i (HfrJ ^ All 

* Javl* (l-®t) (tb— tm) 


(TO) 


and integrate by plotting the inverse cosine function against Au. 

For the total change in lift and moment resulting from this 
correction, equations (56) and (6l (a) and (b)) would be used, with 
the following substitutions; 


for 


Au* and dAu/dtj^ 


for 

a. 

outside the brackets, 0y*/x* 

and ^yb/xb 

for 

a. 

inside the brackets or in Q, 

— 1 and tja 


The total secondary tip correction to the lift (over both halves of 
the wing) would then be (from equation (56)) 
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which should again he integrated by rewriting (BAu/dt^) dt^ as dAu. 
The corresponding correction to the moment can he written similarly 
from equation (6l), The expressions in the brackets become indeter- 
minate when tjn = 0. Their limiting values may be obtained directly 
from the series expansions in Q, since Q(o) = 0. 


Correction for Trailing-Edge Induced Pressures 

The residual pressure at the tip induced by the application of 
the t railing-edge solutions has a somewhat different distribution 
from that induced at the trailing edge by the tip solution. A 
procedure corresponding to that Just described — that is, assuming 
that the pressures Induced by the oblique solutions, as well as by the 
symmetrical solution, all originate at the center of the trailing 
edge — is not so readily Justified as in the preceding case. However, 
the pressures induced by the oblique solutions are quite small, and 
their secondary effect, particularly on the total lift, might in fact 
be neglected, if it were not comparatively simple to take them into 
account with this assumption. 

The approximation is made by substituting in the previously 
given equation (67) for the secondary tralllng-edge correction to the 
local velocity at any point x,y, the value of 


to 


.0.B 

■*fcr«o 


for t a . The approximate induced velocity correction is therefore 

- i _1 cos-i t ftj- t ^ g tptj) (72) 

w Jt 0 =i V-to 

where Au is the total (negative) velocity remaining at the points 
Xjj,s along the tip, found by summing the values from equation (5l) 
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and the integral of equation (50) • 

Equations (29) and (33) are applicable in calculating the 
total correction to the lift and moment, with the following sub- 
stitutions; 

for ua, dAu/dt 0 

for a, except in g(a) and h (a), 0s/xb 

for g(a) and h (a), g(t Q ) and h(t Q ) 

Thus the total secondary trailing— edge change in lift (both halves 
of the wing) would be 

AaL = f ° (a^-Ps) 2 g(t p ) dAu(t p ) (73) 

0&t z t Q =m^ 


and the corresponding correction to the moment can be written 
similarly from equation (33) • 


Further Simplification cf Lift and Moment Calculations 
and Summary of Formulas 


Effect of tip .— The primary corrections to the lift and moment 
due to reducing the pressures at the tip to zero ere obtained without 
difficulty from equations ( 43 .) through ( 47 ) . Calculation of the 
induced pressures at the trailing edge in order to make the secondeiry 
correction is, however, Bomewhat tedious. These pressures vary very 
little over the Mach cone from the tip, and the variation is nearly 
linear. This calculation and, in addition, the integration in 
equation (71) can therefore be eliminated in many cases without any 
significant loss in accuracy by considering Au to be constant at 
some intermediate value between Au* and — u^ (a^. ) , the value at 
the tip. The secondary correction is then simply 

“ ■ ^ (££)"/? (M - « 

where Au may be the average of Au* (equation (30)) and — u^(a^.) 
(equation (10) ) . 

Effect of trailing edge.- Ohe major part of the primary trailing- 


edge correction is accomplished by the application of the symmetrical 
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solution and, in many cases, the integration of equations ( 56 ) and 
( 59 ) can be omitted in calculating the total lift and moment. Then 
the primary effect is given by equation (55). Calculations of the 
secondary effect may be simplified very greatly by ignoring the 
effect of the oblique elements, and using Att as calculated directly 
from equation ( 51 ) in equation ( 73 ) and in the corresponding equation 
for the moment. 


WING WITH CROSS-STREAM TIPS 

The fact, mentioned earlier in this paper, that a drastic 
reduction in lift takes place behind the Mach cone from the tip of 
the highly swept wing suggests immediately the consideration of the 
effect of cutting off the tip areas along a line normal to the stream. 
The resulting plan form (fig. 9) presents obvious structural advan- 
tages and may also be expected to approach the triangular wing in 
its longitudinal stability characteristics. 

In calculating the pressure distribution over such a wing, only 
the trailing-edge solutions apply. For the total lift and moment, 
new integrations must be performed, since the elements of area are 
now different. The characteristic ray a-j- is now replaced by one 
through the inboard end of the cross-stream tip (a* in fig. 9 ) • 

For a ,less than this value. 


d5 ps g 
da 2m? 

and 

X a being defined by equation (48). It is convenient to replace 
Cq by the following expression 


c 0 * 


Then the basic lift is 


( ®t - a i) 
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Lo = Its 2 [ (mt-ax) 2 j 

qa E (y/l-m 2 ) L (mfc^n?) 


L(mt-ai) 


Q-l 2% 


+ COS' 


-1 ai " 

m 


When mt = m. 


— = - *^a- {cos 1 ^ + ~~2 T (2m-ai) \/^-ai 2 - 2 (m-ax) 2 1 1 

la E(/I^) l m 3m 2 L J J 


The decrement of lift due to application of the symmetrical 
solution is (provided, as "before, that the trailing— edge Mach lines 
do not intersect the leading edge) 


at. — 4ai g s 2 

mm^-E (y^l— m? ) 


[§ K (V^l-^t 2 ) - l] 


Except for very high aspect— ratio wings, the trailing edge of 
this configuration will lie entirely in that range over which 
du^/da is negligible compared to u Q . For this reason, the integrar- 
tions necessary to find (AL) a for this case have not been performed. 
They are, however, elementary in nature and can be readily carried 
through by the designer considering a wing of this type. 

The moment of the basic solution about the apex of the wing is 


*k = ^Bs 3 _ Jg 
qa 3mE (v'l-m 2 ) L m t : 


3 (mfc-ai) g j (2ntt z +^) (mt-ai) f Q x g^-aimt 
IJj?) m t 2 ^ n2 \ (mt 2 -*a 2 F^ 2 1_ m(m^-ax) 


jlL ( _J, 

“HiJ \ ^t^ 


mt-ax mf 



+ m ( ttt t~ a i) + m(tnt-ax) 1 + 2 cog -i ax 1 

m-t 2 -^ mt 2 ' J m J 
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When m-fc = m, 


Mo 


S|3s ; 


qa 3 mE (/l— iq 2 ) 

7(ttt-ai) 3 


y^ra 2 — ai 2 1 + 


2 (mr-ai) 

“IS 


1 + i±i) ] 

m / j 


+ cos -1 


15m" 


m 


/ 


(T9) 


The decrement in moment due to the removal of u Q at the 

trailing edge is simply c 0 + — l ^ - Co ) times the decrement 
in lift, or 3 \ m / 

AM . AL 

3inmt 


(80) 


where AL is given by equation (77) * 


Computations have "been made for a ving of this type and will 
he presented in a later section. 


OTHER APPLICATIONS OF THE OBLIQUE SOLUTIONS 


Low Aspect Ratio Wings 

Variations of the method of this paper could he used to obtain 
the pressure distributions for a variety of rectilinear plan forms 
other than the conventional swept wing. For example , the reverse 
delta wing (fig. 10) with trailing edges falling within the Mach 
cones from the tips can be treated by using as the basic solution 
that for a rectangular wing with raked tips (given by Hayes, 
reference 11, and others), the deficiencies from the two— dimensional 
pressure being additive where the Mach cones overlap. The oblique 
solutions of the present report would be applied to cancel the 
pressure deficiency Introduced along the rear of each trailing edge 
by the simple conical field from the opposite tip. Other low-aspect- 
ratio wings can be treated similarly, provided cancellation of pressure 
ahead of a leading edge is not involved. 



h-2 
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Reverse Taper 

Extension of the procedure to wings with reverse taper is 
obvious. In this case the wing with subsonic trailing edge may have 
either subsonic or supersonic leading edge. ('The case of swept— 
forward leading edges requires separate consideration.) Starting 
with the appropriate basic solution, the subsequent procedure and 
formulas of the present report will apply without modification. 


IV - NUMERICAL EXAMPLES 

Hie method has been applied to two wings with 63° sweepback of 
the leading edge, one with considerable taper and one with constant 
chord, at a Mach number of 1.5* The wings are shown in figures 11 
and 12. From the dimensions given, the following quantities, 
required for the computations, can be calculated: 


Parameter 

Tape red 

Untapered 

P 

I.II8 

I.II8 

m 

.570 

.570 

E (y / L^n?) 

I.256 

1.256 

mt 

.906 

.570 

a t 

.528 

• 358 

S 

26 sq ft 

58 sq ft 


It should be noted that the dimensions and calculated values 
omit the tip fairings, which are formed by rotating the (symmetrical) 
tip section about its chord line. The vortex sheet with which a 
lifting surface may be replaced for linear potential flow calculations, 
would have a similar rounding off outboard of the last vortex without 
adding any lift. 

The basic velocity distribution (equation (10)), which depends 
only on m, is the same for both wings. Corrections to this dis- 
tribution, taking place within the Mach cones from the trailing edge 
and tip, cover so small an area of the tapered wing (fig. 11 ), the 
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trailing edge of which almost coincides with the Mach lines, that the 
variation within that area is not of much interest. The effect on 
the total lift and moment will he given later. The pressure distri- 
bution over the constant chord wing will be discussed first in some 
detail. 


PRESSURE DISTRIBUTION OVER SWEPT UNTAPEEED WING 

In practice. It is convenient to compute all velocities in terms 
of u Q , the basic perturbation velocity along the X-axis, in order 
to avoid dime ns ionality and numerous multiplications by a constant 
factor. All but the final figures will therefore show the velocities 
in ratio form. 

The decrement in velocity resulting from removal of constant 
pressure across the wake (sy mme trical solution, equation (51)) is 
plotted against t Q in figure 13- This component is of course 
dependent only on the sweep of the trailing edge relative to the Mach 
lines and would apply equally well to wings of any span and taper 
having the same trailing— edge sweep angle . 

The decrement in velocity due to removal of the remaining 
pressures at the trailing edge was found at various stations span- 
wise and chordwise within the Mach cone from the center of the 
trailing edge and added to the decrement shown for the corresponding 
value of t D in figure 13. Evaluation of this second component 
involved the graphical integration, for each point, of equation (50) 
with respect to a. Figure l4 shows some typical plots of this 
function. The value of a at which each curve falls to zero is a Q 
(equation (53)) for the specified value of x and y. The velocities 
obtained in this way are plotted spanwise for two values of x in 
figure 15. 

The spanwise distribution of velocity decrement within the Mach 
cone from the tip, calculated by equation (37) , is plotted in 
figure l6. 

Figures 13, 15, and 16 represent only the primary corrections. 

In figures 17 (a) to (d) these corrections have been applied to the 
basic velocity distribution along the chord at each of four spanwise 
stations, and t he re sults converted to terms of Ap/qa by the factor 
4u 0 /aV = 4m/fJE (/l-m 2 ) . The magnitude of the extraneous negative 
pressures at the trailing edge remaining after the first approxima- 
tion can be seen in the last two of these figures. Corrections have 
been added to eliminate those errors due to the tip solution by the 
approximate method (equation (70)) of the preceding section. The 
errors introduced by the symmetrical solution have been eliminated 
without resort to approximate methods. The errors induced by the 
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wake gradient solutions were neglected "because of their location 
far "back on the tip and their relatively small values. The effects 
of these secondary corrections are shown in figures 17 (c) and (d). 

For the final curves the pressure distribution has been arbitrarily 
faired to zero, taking cognizance only of the points at which 
succeeding corrections enter. It can be seen by drawing another set 
of reflected Mach lines at the tip of figure 12 that the area affected 
by this approximation is very small and is in that region where viscous 
effects would tend most to invalidate any results calculated by the 
present simplified theory. A three-dimensional view of the final 
pressure distribution is shown in figure 18. 


LIFT AND MOMENT COEFFICIENTS 

The various components of the lift, with the equations from 
which they were obtained, are given in table I for both the tapered 
and constant-chord wing. The last correction was made in two parts, 
one for the symmetrical solution and one for the effect of the 
oblique solutions, assuming the latter to be conical with respect 
to the apex of the trailing edge, in order to estimate the magnitude 
of the error involved in such an assumption. The value given 
represents the upper limit of the possible error. In the secondary 
correction for the tip effect (correction at the trailing edge) 

90 percent of the effect was obtained from the first term of 
equation (71), only 10 percent being contributed by the integral. 

The value of Cj^ obtainable by discarding the trailing-tip 
regions of the constant— chord wing was also estimated, using 
equations (76) and (77)* The modified wing had the same span as the 
original wing, but its area was decreased to 40.8 square feet and its 
aspect ratio increased from I.725 to 2.45 by removal of the trailing 
tips. As a result the value of Cp^ increased from 2.01 to 2.63, 
approaching that for the highly tapered wing. 

(it may be remarked here that the pressure distribution for the 
modified wing is identical with that over the original wing (fig. 17) 
up to the point at which the tips are cut off, where the pressure 
drops discontinuously to zero.) 

The components of pitching moment for the three wings are 
presented in table II. Second— approximation corrections were not 
made where the corresponding component of lift was negligible. 


DRAG DUE TO LIFT 

With the value of calculated, it is possible to find 

immediately the drag due to lift for the specified wings. As pointed 
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out in reference 12, where similar calculations are made for a more 
limited class of wings, the flat wing at an angle of attack is 
subject to a drag which is merely the component of the lift in the 
flight direction, or a times the lift, and to a thrust which is 
due to the infinite suction force at the leading edge. If, as in 
the examples calculated, the Mach lines from the trailing edge do 
not cross the leading edge, the thrust on the leading edge is 
exactly that on the triangular wing of equal span at the same angle 
of attack. From reference 12, this thrust, in coefficient form, is 

py/l-m 2 2 

4itm tL A 


where Cl^ is the lift coefficient of the triangular* wing. From 
reference 6, of the triangular wing with m<l is 

Using this value and the area of the triangular wing ps 2 /m, 

pE (\Zl-m 2 ) 

the thrust force can he expressed in terms of the sweep of the trap- 
ezoidal wing as follows: 


T = ffaVl^g qaJ 2 
[E 4/l-m?) ] 2 


The total drag for the swept wing is therefore 
qa 2 Cr S - T = qa 2 \ C t S 

^ L ^ [E(v ^^)] 2 


and the drag coefficient may he written 

Cd = icr (&) [a 2 

[E(v^^)] 2 \S/J 

or 

Qd _ 1 i c _ ( a£ \ l 

CiF " CL a 2 1 ^ [E (v/l-m 2 )] 2 \ S / J 


(81) 


(82) 


For the three wings described, the results may he summarized as 
follows : 
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0.160 
.323 

.253 

CONCLUDING REMARKS 

The calculation of the theoretical load over a swept-back wing 
of moderate aspect ratio at a Mach n umb er not too close to one 
appears to he feasible by the method described in this paper. 
total lift can be approximate d within 5 percent for a considerable 
range of practical plan forms with only the first step of the 
successive approximations outlined. In the cases calculated, the 
major effects were obtained with only one graphical integration for 
each wing (that for the primary tip effect). The moment is not as 
readily approximated; son© secondary corrections will generally prove 
necessary. 

Although in the examples calculated the oblique solutions con- 
tributed only about 1 percent to the final lift and 2 percent to the 
moment, they are expected to prove useful for other cases, particu- 
larly for low- aspect -ratio wings of other than rectangular plan form. 

At Mach numbers close to 1.0, the method is no longer practical 
especially for predicting pressure distributions. It appears that 
some other line of approach is required for this case. 

The results of the calculations emphasize the indications of the 
theory that the determining factor in the aerodynamic performance of 
wings, as far as differences in plan form are concerned, is the dis- 
position of area relative to the Mach lines. The apparent beneficial 
effect of taper and increased aspect ratio, for example, is largely 
the result of reducing the percentage of area within the Mach cones 
from the tip and trailing edge. The negative lift calculated for the 
tip areas, in particular, may be expected to have a serious effect 
on the longitudinal stability variation with Mach number. With the 
proper thickness distribution, the wing with cross-stream tips may 
provide the solution to this problem. 
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Wing 

Aspect 

ratio 

_a 

Tapered wing 

3.85 

3.09 

Untapered wing 

1.73 

2.01 

Modified untapered 



wing 

2.45 

2.63 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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APPENDIX A 
Summary of Symbols 

x,y,z cartesian coordinates in the stream direction, across the 
stream, and in the vertical direction, respectively 

V free- stream velocity 

u,v,v perturbation velocities in the stream direction, across the 
stream, and in the vertical direction, respectively 

M free— stream Mach number 


V^-l 

density of air 


dynamic pressure 


pressure difference between upper and lower surfaces, or 
local lift 

root chord 




semi span 
wing area 

angle of attack, radians 

angle of sweep of the leading edge 

slope of leading edge _ ^ cQt A 
slope of Mach lines 

slope of trailing edge 
slope of Mach lines 

slope of any ray through origin _ 0y 
slope of Mach lines ~ x 

Also used to designate a ray through the origin or a 
constant— load element having such a ray as one side, the 
other side being a boundary of the wing surface. Used 
as a subscript to designate any quantity pertaining to 
such an element. 
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a t 

x a#7a 

x Td^ 

x* )7 * 

t a 

% 
to 
tm 
t* ' 

m 

Hi 


value of a corresponding to most rearward element affecting 
a specified point x,y. On tip, a D = -- 


x+p (y— s ) 


On trailing edge, a Q = mfc 


py+c 0 -x 

Py+Comt-x 


slope of ray through traillng-edge tip _ ps 

slope of Mach lines Co+Ps/mt 

coordinates of apex of any constant-load element a 

coordinates of apex of elements used in secondary corrections 

coordinates of the point at which the Mach line from the tip 
crosses the trailing edge 

slope of ray through apex of element a _ ^ y— y ft a 

slope of Mach line. x— x a 

point x,y 

slope of ray through Xk,yh 


= 3 


y-7h 


slope of Mach lines x-3C b 

slope of ray through trailing— edge apex 

slope of Mach lines 

> 

slope of ray from leading-edge tip y_g 

slope of Mach line 

slope of ray through x*,y* _ ^ y— y* 
slope of Mach lines x— x* 


x— c f 


- M 
m 


hasic perturbation velocity as given hy solution for 
triangular wing 

value of ua for wing with supersonic leading edge, 
when a > 1 


u D value of ua for wing with subsonic leading edge when 

a=0 (center line of wing) 

u a incremental velocity on any constant— load element a 

Au correction to the basic incremental velocity at any point 

(Au) a correction to the basic incremental velocity induced by a 

single element a 



MCA TN Wo. 1555 


k9 


Au* Decrement in velocity u induced at x*,y* by cancellation 
of pressures at the tip 

L q basic lift for entire wing 

Mo basic moment 

(£L) a correction to total lift due to application of one constant- 
load element 


AL 
<D ‘ 

9 

T 

U 

u 

c 


c 

g(a) 

h(a) 

Q 


correction to the basic lift 


complex velocity potential function 
real part of 0 
imaginary part of 4 

complex velocity function of which u is the real part 


imaginary part of U 

, argunent for solutions of the supersonic 

flow equation 


y+iz 

/x 2 -P 2 y a -0z 2 




1 / a+a 2 

2 (l+m t ) (mf-a) \j mt+mt® 


a(l+mt) 



E elliptic integral of the second kind 

F elliptic integral of the first kind 

K complete elliptic integral of the first kind 
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APPENDIX B 

Integration for Loss of Lift at the Tip 
of the Highly Swept Wing 

From equations (42), (44), and (10) 


where 


O’ (a) 


AL _ 4mt 2 m 2 s 2 P m 

qa £ (/l-m 2 ) J& £ 


l/m^-a 2 


da 


a-at 

at 2 a 2 (mt-a) 


( 


a t . m f~ a t \ ( I a-fa 2 a_ \ 

a mf-a J \ ymfc+mt 2 mt J 


(Al) 


ar-at 

2^111^ +m^ 5 ~ /a+a 2 


(A2) 


The terms in O’ (a) are of two types; namely, those that 
contain i/ a+a g a nd those that do not. The former combine with the 
radical i/mP-a 2 in (Al), to form elliptic integrals of the first, 
second, and third kinds. The latter give rise to terms in 
equation (Al) which are immediately integrahle by elementary means. 
It is convenient, therefore, to consider the integral in two parts, 
writing 



&» (a) 

^/m 2 — a 2 


da = Ii + I 2 


where Ii is that part of the integral not requiring elliptic 
integrals . 


Then 


rm 

\"L 


t 


. (*± + 2££k \ 

, t at 2 mta(mt-a) \a n%-a J 


da 


3 f ooa _ 1 2tairjf_ + gf^tg- t ^5=^i 

a t 2jn t (mfcS-m 2 ) L y' mt 2_ m 2 m(m-t-a-t) 


(A3) 
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Die remaining terms, involving /a-m 2 and y/m 2 — a 2 , are inte- 
grated by means of the substitution 3 



so that 



The terms involving require in addition the following func- 

tions and parameters: 


n 


m ( m f*~l) 

mt-m 


and a = tn 1 



(A5) 


In the final formula for I2, given below, the elliptic 
functions have been re-expressed, as far as possible, in terms of 
the original parameters, in order to simplify the computing pro- 
cedure* The result is 


The symbols u and a used throughout the remainder of this 
appendix are the standard notation for elliptic functions and 
are not to be confused with the aerodynamic symbols of the text. 
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T- - 1 ( _^h 1 ] F( , 

mt/anm,. (l+m^) \^L \ mfc+1 m/ m 2 J 


+ f 1 ££££ 'Wq^k) + ( L jg-pE 

m L \ “t / V 6 11 % / V 2na t (at+l) 


(14m) (mt-eut ) 8 


, du m 

1 + r-S- + — -*- 


1 -HDtt mt 


L 

-Hat). 


n (n,k,<p) 


+ 2m(i4ffl) dilate, k,q>) 
mf-m Si ~ 


cp = sin' 


1 p-at 


i(a t +l) 


n - r K> — 4 a 

3 in l+n sn g u 


Is Legendre’s form of the elliptic integral of the third kind. 
In the present problem the integral 


■ui (at) 


du 

(l+n sn 2 u ) 2 


alBO arises. Differentiation of II 3 with respect to the 
parameter n gives 


^_ 3 =_ r l 

dn Jr, 


sn?u du 
(l+n sn 2 u ) 2 


.1/*' 

n Jo 


du 

(l+n sn 2 u ) 2 


-1 / 

n Jo 


u(at) 


l+n sn 2 u 


which may then he solved for 


du 

(l+n sr^u ) 2 
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This relation is the source of the term in Sn 3 /Sn in equation (a 6 ). 
While II 3 and its derivative are not tabulated, they can he eval- 
uated for individual cases by the following method: 


We first note that n>o. The procedure to evaluate II 3 
required in that case has been outlined in some detail as case XV 
in reference 13 . The following is an adaptation of that procedure 
to make use of available tabulated functions. The parameter a 
(equation (A 5 ) ) is introduced. Then, if u-^ is the value of u 
when a = a^. 


n 3 (n,k,cp) = ut 



sn (a,k* )cn (a,k* ) 

dn(a,k’) 


iZ (ia) + ~ 
u t 



with 

iZ (ia) = -2SL. + E (a,k») - a H fa,*' 

2 KK* K‘ cn (a,k* ) 


and 


^ = tan 


_t 2q sin2a n ut sinh 2a n q — 2q 4 sinka n ut ainh4 a n a+... 

1— 2q cos2a n ut cosh 2a n a + 2q 4 cos4a n ut cosh 4a n ct-. . . 


The elliptic functions sn(a,k’), cn(a,k*), dn(a,k*) can be computed 
from the value of tn(a,k* ) (equation (A5)), while u^ is simply 
F (T,k) . In Z(ia), K,K f and E,E’ are the complete elliptic 
integrals of moduli k and k*, respectively. In ^ , q is the 

— itK» 

El 

quantity e , and a n is the asymptotic value of the sequence 
defined by the recurrence formulae a n = i- + bj^) and 

tn = \/ a nr-i tn-x with ao = 1 and b Q = k* . (See reference 13 . ) 


The derivative of II 3 with respect to n is 


1 3 WAUU x oapou u uu u 18 Which 

may be obtained for this case (n > o) in the form 

fa.it') ( rae. fo.fc . *) _ an g (a,k')l { r.fek.) 

Sn 2n dn(a,k* ) \ Ldn 2 (a,k*) J l L 

\ 



u + \j/ 



sn (a,k* )cn (a.k* ) 
dn(a,k') 




u 
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Substitution of the limits and the expressions in equation (A4), 
gives 


5n 3 (n,9,k) _ Bn (a,k* )cn (a,k f ) 
dn 2n dn (a,k* ) 



sn (a,k* )cn (a,k* ) 
dn(a,k*) 


]{[*«■*■> 

+ 

| j^ l + dn 2 (a,k* ) J F(<P,k) 


- 1 (9,k) 


m-t+1 


/ a t (m 2 — a-^ 2 ) 
V 2m(a t +l) 
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V////M Region of 

constant 

pressure 



Figure ir Conical flow regions on wings with various 
degrees of sweep relative to the Mach lines. 



Fig. 2 
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Figure 2r Detail of tip of tow aspect ratio wing with 
supersonic leading edge , showing axes and 
symbols used in formulas ,and basic velocity distribution. 





Fig. 4 
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Fig. 5 



Figure 5r Detail of trailing edge of highly swept wing showing 
oblique constant -lift element and induced velocity 
distribution. 



Fig. 6 
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Figure 6.- Regions of tip and trailing edge interaction. 
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Fig. 7 




Figure 7 -Second step of successive approximations. Correction 
for pressures at tip induced by traiiing-edge solution. 
Shaded area denotes extraneous pressure field. 




Fig. 8 
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Figure dr Approximate cancellation of extraneous pressures 
behind trailing edge by single conical field from the 
leading edge Up. 




Fig. 10 
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Figure iO.— Reverse delta wing . 
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Fig. 15 



Figure 15- Two typical spanwise incremental velocity distributions, representing 
modifications to the basic solution within the Mach cone from the 
trailing edge apex (contribution of oblique solutions only). 
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Fig. 16 



Figure !6r~ Spanwise distribution of incremental 
velocity at various chordwise stations 
within the Mach cone from the tip of the 
untapered wing. 




(b) £ =.5/3 
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Figure 17 r Chord wise pressure distributions calculated 
for the untapered wing. 
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Figure !7- (continued) 
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Figure /7.— (concluded) 







